LATTICE-BASED UNSUPERVISED MAXIMUM LIKELfflOOD LINEAR 
REGRESSION FOR SPEAKER ADAPTATION 

Field of the Invention 

The present invention generally relates to methods and arrangements for providing 
speaker adaptation in connection with speech recognition. 

Background of the Invention 

Acoustic adaptation is playing an increasingly important role in speech recognition 
systems, to compensate for the acoustic mismatch between training and test data, and also 
to adapt speaker-independent systems to individual speakers. Most speech recognition 
systems use acoustic models that include multi-dimensional gaussians that model the 
probability density function (pdf) of the feature vectors for different classes. (For general 
background on speech recognition, including gaussian mixture pdf s, see, e.g. 
Fundamentals of Speech Recognition^ Lawrence Rabiner and Biing-Hwang Juang, 
Prentice Hall, 1993; and Statistical Methods for Speech Recognition, Frederick Jelinek, 
The MIT Press, 1997.) A commonly used adaptation technique in this connection is 
maximum likeUhood linear regression (MLLR), which assumes that the parameters of the 
gaussians are transformed by an affine transform into parameters that better match the test 
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or adaptation data. In a simple implementation, the mean Ui of each gaussian gi is 
transformed according to Ui' = Aui where A is the transform matrix, and ui is optionally 
padded with ones to represent an offset. The transform is chosen so as to maximize the 
probability of a collection of adaptation data with associated transcriptions. In more 
sophisticated implementations, the gaussian variances may also be adjusted. MLLR is 
further discussed, for instance, in Leggetter et al., "Speaker Adaptation of Continuous 
Density HMM's Using Multivariate Linear Regression", Proceedings of ICSLP '94, 
Yokohama, Japan, 1994. This technique is also often used in "unsupervised" mode, where 
the correct transcription of the adaptation data is not known, and a first pass decoding 
using a speaker independent system is used to produce an initial transcription. 

Although MLLR appears to work fairly well even when the unsupervised 
transcription is mildly erroneous, it is recognized herein that fiirther improvements are 
possible. 

Accordingly, a need has been recognized, inter alia^ in connection with improving 
upon the shortcomings and disadvantages associated with conventional arrangements such 
as those discussed above. 
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Summary of the Invention 



In accordance with at least one presently preferred embodiment of the present 
invention, it is presently recognized that it is possible to improve upon the performance of 
MLLR, even when the unsupervised transcription is mildly erroneous, by taking into 
account the fact that the initial transcription contains errors. This may be accomplished, 
for example, by considering not just the "1-best" (i.e., single best) transcription produced 
during the first pass decoding, but the top N candidates. (See, for example, Jelinek 
[1997], supra, for a description of "N-best" decoding and definitions associated 
therewith.) Alternatively, if the first pass decoding produces a word graph, this can be 
used as the reference word graph, instead of the 1-best or N-best reference transcriptions. 
In contrast to an N-best list, which simply enumerates a relatively small number (e.g. 100 
or 1000) of likely word sequences, a word graph is a compact representation of all the 
word sequences that have any appreciable probability. An example of a word graph is 
illustrated in Figure 2. 

Broadly contemplated herein, in accordance with at least one presently preferred 
embodiment of the present invention, is a formulation that affinely transforms the means of 
the gaussians to maximize the log likelihood of the adaptation data under the assumption 
that a word graph is available that represents all possible word sequences that correspond 
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to the adaptation data. The word graph is produced during a first pass decoding with 
speaker independent models. It is also possible to consider only those regions of the word 
graph that represent a high confidence of being correct to further improve the 
performance. 

In one aspect, the present invention provides a method of providing speaker 
adaptation in speech recognition, the method comprising the steps of: providing at least 
one speech recognition model; accepting speaker data; generating a word lattice based on 
the speaker data; and adapting at least one of the speaker data and the at least one speech 
recognition model in a manner to maximize the likelihood of the speaker data with respect 
to the generated word lattice. 

In another aspect, the present invention provides an apparatus for providing 
speaker adaptation in speech recognition, the apparatus comprising: at least one speech 
recognition model; an accepting arrangement which accepts speaker data; a lattice 
generator which generates a word lattice based on the speaker data; and a processing 
arrangement which adapts at least one of the speaker data and the at least one speech 
recognition model in a manner to maximize the likelihood of the speaker data with respect 
to the generated word lattice. 
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Furthermore, in another aspect, the present invention provides a program storage 
device readable by machine, tangibly embodying a program of instructions executable by 
the machine to perform method steps for providing speaker adaptation in speech 
recognition, the method comprising the steps of providing at least one speech recognition 
model; accepting speaker data; generating a word lattice based on the speaker data; and 
adapting at least one of the speaker data and the at least one speech recognition model in a 
manner to maximize the likelihood of the speaker data with respect to the generated word 
lattice. 

For a better understanding of the present invention, together with other and further 
features and advantages thereof, reference is made to the following description, taken in 
conjunction with the accompanying drawings, and the scope of the invention will be 
pointed out in the appended claims. 

Brief Description of the Drawines 

Fig. 1 illustrates a Hidden Markov Model (HMM) structure used to generate 
Maximum A-Posteriori Probability (MAP) lattices; 

Fig. 2 illustrates word traces produced by the MAP lattice HMM, and their 
connection into a word lattice; 
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Fig. 3 illustrates a histogram of state posterior probabilities; and 

Fig. 4 illustrates a graph of word error rate versus confidence threshold. 

Description of the Preferred Embodiments 

Throughout the present disclosure, various terms are utilized that are generally 
well-known to those of ordinary skill in the art. For a more in-depth definition of such 
terms, any of several sources may be relied upon, including Rabiner and Juang (1993), 
supra, and Jelinek (1997), supra. First, a theoretical framework is discussed in connection 
with at least one presently preferred embodiment of the present invention. 

Also, the article "Lattice Based Unsupervised MLLR for Speaker Adaptation in 
Speech Recognition Systems" (Mukund Padmanabhan, George Saon, Geoffrey Zweig, 
ISCA ITRW ASR2000 Paris, France 2000 [http://]www-tlp,limsi.fi-/asr2000) is hereby 
fiiUy incorporated by reference as if set forth in its entirety herein. 

In a typical speech recognition system, the speech signal is represented as a 
sequence of observation vectors, and throughout the present discussion, yt denotes the 
multi-dimensional observation at time andy;^ denotes the 7 observations corresponding 
to the adaptation data. The pdf s of each context dependent phonetic state s is modeled 
by a single gaussian (this can be easily generalized to mixtures of gaussians) with mean 
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and diagonal covariance iLis , As, 9 \s used to indicated the current values of the gaussian 
parameters, and ^is used to denote the future (adapted) values to be estimated. The 
probability density of the observation j^r given the pdf of state s is denoted p$(yt/s). It 
will presently be assumed that ^and 6 are related in the following way: = Aju^, = 

Ay, i.e., only the current means of the gaussians are linearly transformed, and all means are 
transformed by the same matrix A. 

Typically, in an MLLR framework, the general objective is defined as follows: 
given a transcription w of the adaptation data, find 9 (or equivalently A) so that the log 
hkelihood of the adaptation data, y/ is maximized. The transcription w can be represented 
as a sequence of K states 5; . Sk, and the T observation frames can be aligned with this 
sequence of states. However, the alignment of the T frames with the sequence of states is 
not known. Let St denote the state at time The objective is to find the maximum 
likelihood transform theta, and can now be written as: 

^ * = arg max log \ (/ )1 

s ^ 

= argmax ^ log [/j^- (1) 
= argmax Y,Pe (^1 )log[/'e (•>'^^0] 
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In a lattice-based MLLR currently contemplated in accordance with at least one 
embodiment of the present invention, it is assumed that the word sequence, and thus the 
state sequence si^ , corresponding to the adaptation data cannot be uniquely identified, 
and this uncertainty is incorporated in the form of a lattice or word graph. Preferably, the 
word graph is produced by a first pass decoding with speaker independent models. The 
formulation of the maximum likelihood problem is essentially identical to equation (1), but 
with one significant difference. In (1), the states St were assumed to belong to the alphabet 
of K states Si . . . Sk, with the only allowed transitions being Si and Si -> Sf+L In a 
lattice-based MLLR formulation according to at least one presently preferred embodiment 
of the present invention, the transition between the states is dictated by the structure of the 
word graph. Additionally, it is possible to take into account the language model 
probabilities (which are ignored in the MLLR formulation), by incorporating them into the 
transition probability corresponding to the transition from the final state of a word in the 
word graph to the initial state of the next connected word in the word graph. 

The disclosure now turns to a decoding strategy for producing word graphs. 

In accordance with at least one presently preferred embodiment of the present 
invention, a Maximum A-Posteriori Probability (MAP) word lattice is preferably generated 
using word internal acoustic models and a bigram language model. MAP lattices and 
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bigram language models are discussed generally in several publications, including (Jelinek, 
1997). To construct the lattice, it may be assumed that the utterance in question is 
produced by an HMM with a structure such as that shown in Figure 1 . Each 
pronunciation variant in the vocabulary appears as a linear sequence of phones in the 
HMM, and the structure of this model permits the use of word-internal context dependent 
phones. Preferably, a bigram language model is used with modified Kneser-Ney 
smoothing {see, for example, Kneser and Ney, "Improved Backing-oflF for n-gram 
Language Modeling", Proceedings of the International Conference on Acoustics, Speech, 
and Signal Processing, 1995) . Here, there is an arc fi-om the end of each word to a null 
word-boundary state, and this arc has a transition probability equal to the back-off 
probability for the word. From the word-boundary state, there is an arc to the beginning 
of each word, labeled with the unigram probability. For word pairs for which there is a 
direct bigram probability, an arc is preferably introduced fi-om the end of the first word to 
the beginning of the second, and this arc preferably has a transition probability equal to the 
discounted bigram probability. Preferably, the dynamic range of the acoustic and 
language-model probabilities is normalized by using an appropriate language model 
weight, such as 15. 

Preferably, the MAP lattice is constructed by computing the posterior state 
occupancy probabilities for each state at each time: 
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where al = P(y}\ St=s) and fil = Pfy^/ I St=s) , and then computation posterior word 

occupancy probabilities by summing over all the states interior to each word. That is, if 
w. is the set of states in word Wi , then the following is preferably computed at each time 

frame: 



Preferably, the likeliest words are kept track of at each frame, and these are preferably 
output as a first step in the processing. 

It will be noticed that a word will be on the list of "likeliest words" for a period of 
10 time, and thereafter will "fall off" that list. Thus, the output of the first step may 

preferably be a set of word traces, as illustrated in Figure 2. The horizontal axis is time, 
while the vertical axis ranges over all the pronunciation variants. 

Preferably, the next step will be to connect the word traces into a lattice. Many 
connection schemes are possible, but it has been found that the following strategy is quite 
15 effective. It requires that one more quantity be computed as the word traces are 
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generated: the temporal midpoint of each trace as computed from the first moment of its 
posterior probability: 



i^t- start * V / 

St~end ^ 
t-start 



To construct an actual lattice, a connection is preferably added from the end of one 
5 word trace to the beginning of another if the two overlap, and the midpoint of the second 

s 

H is to the right of the midpoint of the first This is illustrated at the bottom of Figure 2. (It 

P 

J 2 has also been found to be convenient to discard traces that do not persist for a minimum 

period of time, or which do not reach an absolute threshold in posterior probability.) 

ry To evaluate the lattices, the oracle worderror rate (i.e. the error rate of the single 

s 

O 10 path through the lattice that has the smallest edit distance from the reference script) can be 

JSCS; 

computed. This is the best worderror rate that can be achieved by any subsequent 
processing to extract a single path from the lattice. For voicemail transcription, the MAP 
lattices have an oracle word error rate of about 9%, and the ratio of the number of word 
occurrences in the lattices to the number of words in the reference scripts is about 64. 
15 Due to the rather lax requirements for adding links between words, the average indegree 
for a word is 74; that is, there are about 74 possible predecessors for each word in the 
graph. The MAP lattice that is produced in this way is suitable for a bigram language 
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model: the arcs between wordends can be labeled with bigram transition probabilities, but 
is too large for a straightforward expansion to trigram context. In order to reduce its size, 
a second pass is preferably made, where the posterior probability of transitioning along the 
arcs that connect wordtraces is computed. That is, if is the last state in one word trace 
and Sj is the first state in a successor and is the weighted language model transition 
probability of seeing the two words in succession, one may compute: 



The above equation represents the posterior probability of being in state 5/ at time t 
and in state Sj at time / + 7, and transitioning between the words at an intermediate time. 
For each link between word traces, this quantity is summed over all time to get the total 
probability that the two words occurred sequentially; the links with the lowest posteriors 
are then discarded. It should be noted that a separate quantity is preferably computed for 
every link in the lattice. Thus, even if two links connect traces with the same word labels, 
the links will in general receive different posterior probabilities because the traces will lie 
in different parts of the lattice, and therefore tend to align to different segments of the 
acoustic data. 
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As in Mangu et al., "Lattice Compression in the Consensual Post-Processing 
Framework" (Proceedings of SCI/IS AS, Orlando, Fla., 1999), it has been found that over 
95% of the links can be removed without a major loss of accuracy. Here, it was found 
that pruned lattices had an average indegree a little under 4, and an oracle error rate of 
about 11%. After pruning, lattices were expanded to a trigram context, and the posterior 
state occupancy probabilities needed for MLLR were computed with a modified Kneser- 
Ney trigram language model, along with leftword context dependent acoustic models. 

Accordingly, the disclosure now turns to a discussion of a confidence-related 
pruning method that enables regions of low confidence to be discarded. 

Word lattices have been used in a variety of confidence estimation schemes (see, 
for example, Kemp et. al., "Estimating Confidence Using Word Lattices" (Proceedings of 
ICASSP '97, 1997) and Evermann et al, "Large Vocabulary Decoding and Confidence 
Estimation Using Word Posterior Probabilities (Proceedings of ICASSP '00, 2000). Here, 
the simplest possible measure posterior phone probability was explored for discarding 
interpretations in which there was low confidence. It is to be recalled that, as a first step 
in MLLR, the posterior gaussian probabilities are computed for all the gaussians in the 
system. This is computed on a phone-by-phone basis, first computing the posterior phone 
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probability, and then multiplying by the relative activations for the gaussians associated 
with the phone. For phone Si with gaussian mixture , and for a specific time frame yt , 



Since the gaussian posteriors are used to define a set of Unear equations that are 
solved for the MLLR transform, it is reasonable to assume that noisy or uncertain 
estimates of the posteriors will lead to a poor estimate of the MLLR transform. To 
examine the truth of this hypothesis, the MLLR transform was estimated from subsets of 
the data, using only those estimates of P(St ^ Si\y/) that were above a threshold, typically 
0.7 to 0.9. 

The experiments were performed on a voicemail transcription task. (For a general 
discussion of voicemail transcription, see Padmanabhan et al., "Recent Improvements in 
Voicemail Transcription" (Proceedings of EUROSPEECH '99, Budapest, Hungary, 1999). 
The speaker independent system has 23 13 context dependent phones, and 134,000 
diagonal gaussian mixture components, and was trained on approximately 70 hours of 
data. The feature vectors are obtained in the following way: 24 dimensional cepstral 
vectors are computed every 10ms (with a window size of 25ms). Every 9 consecutive 
cepstral vectors are spliced together forming a 216 dimensional vector which is then 
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projected down to 39 dimensions using heteroscedastic discriminant analysis and 
maximum likelihood linear transforms (see Saon et al., "Maximum Likelihood 
Discriminant Feature Spaces", to appear in Proceedings of ICASSP '2000, Istanbul, 2000). 

The test set contains 86 randomly selected voicemail messages (approximately 
7000 words). For every test message, a firstpass speaker independent decoding produced 
a MAP word lattice described in section 3. For the MLLR statistics we used phone and 
gaussian posteriors as described in section 4. The regression classes for MLLR were 
defined in the following way: first all the mixture components within a phone were 
bottom-up clustered using a minimum likelihood distance and next, the representatives for 
all the phones were clustered again until reaching one root node. The number of MLLR 
transforms that will be computed depends on the number of counts that particular nodes in 
the regression tree get. In practice, a minimum threshold of 1500 was found to be useful. 
For voicemail messages which are typically 10 to 50 seconds long this results in 
computing 13 transforms per message. 

Figure 3 shows the histogram of the non zero phone posteriors computed over all 
the test sentences. It is to be noted that, first, there are a significant number of entries with 
moderate (0.1 - 0.9) probabilities. Secondly, although there are a significant number of 
entries at the leftend of the histogram, they have such low probabilities that they account 
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for an insignificant amount of probability mass. This suggests that one can use high values 
for the confidence thresholds on the posteriors without losing too much adaptation data. 

Figure 4 shows the word error rate as a function of the confidence threshold. The 
optimal results were obtained for a threshold of 0.8. Increasing the threshold above this 
value results in discarding too much adaptation data which counters the effect of using 
only alignments in which one is very confident. 

Finally, Table I (herebelow) compares the word error rates of the speaker 
independent system, 1-best MLLR, lattice MLLR and confidence-based lattice MLLR. 
The overall improvement of the confidence-based lattice MLLR over the Ibest MLLR is 
about 1.8% relative and has been found to be consistent across different test sets. It is 
expected that the application of iterative MLLR, i.e. repeated data alignment and 
transform estimation, will increase the improvement. This is because the lattice has more 
correct words to align to than the Ibest transcription. For comparison, Wallhoff et al,, 
"Frame-Discriminative and Confidence-Driven Adaptation for LVCSR" (Proceedings of 
ICASSP '00, 2000) cites a gain on a "Wall Street Journal" task of 34% relative over 
standard MLLR by combining confidence measures with MLLR. 
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TABLE 1 



System Word Error Rate 

Baseline (SI) 33.72% 

1-bestMLLR 32.14% 

Lattice MLLR 3L98% 

Lattice MLLR + threshold 3 1 .56% 



In recapitulation, the present invention, in accordance with at least one presently 
preferred embodiment, broadly contemplates the use of a word lattice in conjunction with 
MLLR. Rather than adjusting the gaussian means to maximize the likelihood of the data 
given a single decoded script, a transform was generated that maximized the likelihood of 
the data given a set of word hypotheses concisely represented in a word lattice. It was 
found that the use of a lattice alone produces an improvement, and also that one can gain a 
more significant improvement by discarding statistics in which one has low confidence. 
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In further recapitulation, it will be appreciated from the foregoing that the use of 
lattice-based information for unsupervised speaker adaptation is explored herein. It is 
recognized that, as initially formulated, MLLR aims to linearly transform the means of the 
gaussian models in order to maximize the likelihood of the adaptation data given the 
correct hypothesis (supervised MLLR) or the decoded hypothesis (unsupervised MLLR). 
For the latter, if the first-pass decoded hypothesis is significantly erroneous (as is usually 
the case for large vocabulary telephony applications), MLLR will ofl:en find a transform 
that increases the likelihood for the incorrect models, and may even lower the likelihood of 
the correct hypothesis. Since the oracle word error rate of a lattice is much lower than 
that of the 1-best or N-best hypothesis, by performing adaptation against a word lattice, 
correct models are more likely to be used in estimating a transform. Further, a particular 
type of lattice proposed herein enables the use of a natural confidence measure given by 
the posterior occupancy probability of a state, that is, the statistics of a particular state will 
be updated with the current frame only if the a posteriori probability of the state at that 
particular time is greater than a predetermined threshold. Experiments performed on a 
voicemail speech recognition task indicate a relative 2% improvement in the word error 
rate of lattice MLLR over 1-best MLLR. 
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The present invention is applicable to all particular forms of MLLR, including 
those in which the gaussian variances are transformed, and those in which the feature 
vectors are transformed. 

It is to be understood that the present invention, in accordance with at least one 
presently preferred embodiment, includes at least one speech recognition model, an 
accepting arrangement which accepts speaker data, a lattice generator which generates a 
word lattice based on the speaker data, and a processing arrangement which adapts at 
least one of the speaker data and the at least one speech recognition model in a manner to 
maximize the likelihood of the speaker data with respect to the generated word lattice. 
Together, the accepting arrangement, lattice generator and processing arrangement may 
be implemented on at least one general-purpose computer running suitable software 
programs. These may also be implemented on at least one Integrated Circuit or part of at 
least one Integrated Circuit, Thus, it is to be understood that the invention may be 
implemented in hardware, software, or a combination of both. 

If not otherwise stated herein , it is to be assumed that all patents, patent 
applications, patent publications and other publications (including web-based publications) 
mentioned and cited herein are hereby fixUy incorporated by reference herein as if set forth 
in their entirety herein. 
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Although illustrative embodiments of the present invention have been described 
herein with reference to the accompanying drawings, it is to be understood that the 
invention is not limited to those precise embodiments, and that various other changes and 
modifications may be affected therein by one skilled in the art without departing from the 
scope or spirit of the invention. 
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